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Terahertz  Wave  Interaction  with  Biological  Macromolecules 
(Experiment  and  Simulation) 


Ohjective: 

To  discover  and  understand  the  fundamental  physical  principles  governing  the 
interactions  of  biological  molecules  with  eleetromagnetie  radiation  in  the  terahertz  region 
of  the  speetrum.  The  primary  goals  are  to  observe  the  resonant  modes  in  absorption  spectra 
within  the  submillimeter  wave  range  and  to  demonstrate  the  fundamental  character  of  the 
submillimeter-wave  spectral  absorption,  thereby,  establishing  a  physical  confirmation  that  very 
long  wavelength  resonant  features  result  from  electromagnetic  wave  interactions  with  the  DNA 
and  RNA  macromolecules  -  phonon  modes. 

Impact  of  the  Research: 

The  lead  research  component  is  focused  on  the  molecular  structure,  dynamical  properties,  and  the 
THz  interaction  mechanisms  related  to  CBW  agents  in  order  to  develop  a  scientific  base  for  THz- 
frequency  sensing.  The  impact  of  this  research  will  be  to  establish  the  initial  foundation  for  the 
future  use  of  terahertz  spectroscopy  in  the  identification  and  characterization  of  biological 
macromolecules  and  related  materials. 

Technology  Transfers  and  Research  Collaborations: 

This  research  results  from  collaboration  with  Dr.  D.  Woolard  (ARO),  Dr.  A.  Samuels  and  Dr.  J. 
Jensen  (Edgewood  Chemical  and  Biological  Center),  with  the  group  of  Prof  T.  Crowe  (UVA  and 
Virginia  Diodes  Inc.),  groups  of  Prof  R.  Weikle  (UVA),  Prof.  Lucas  Tamm  (UVA),  Prof  E. 
Fernandez  (UVA),  Prof.  D.  Theodoresku  and  Prof  H.  Frierson  (UVA),  the  group  of  Prof  N. 
Stewart  (the  University  of  Tennessee),  with  the  group  of  Prof  E.  Brown  (UCLA),  Lowell 
University  in  Massachusetts,  Goodrich  Co.  Results  of  characterization  immediately  transfer  to 
Edgewood  Chemical  and  Biological  Center.  Developed  technology  for  fingerprinting  of 
biological  materials  in  the  THz  spectral  range  was  used  to  construct  a  predicted  database  for 
THz  signature  of  bio-warfare  agents  using  experimental  transmission  speetroscopy  and 
modeling  in  a  DTRA/CBT  Project  ‘Wew  Concepts  for  Detection  Biological  Targets: 
Terahertz  Signature  Data  Base  Generation“,the  proposal  number;  P-472 13_EL-CDP, 
grant  No:  W91  lNF-05-1-0033  (BAA05DET014:  UVA,  Task:  47213-EL-CDP,  Dr.  Ngai 
M.  Wong  ),  PI  T.  Globus.  The  results  will  be  also  used  in  future  development  of  biosensing 
technology  (Virginia  Diodes). 

Summary  of  Achievements: 

•  The  fundamental  physical  nature  of  the  observed  resonance  structures  in  THz  spectra  of 
biological  macromolecules  which  is  caused  by  the  internal  vibration  modes  in 
macromolecules  was  further  confirmed  and  studied.  The  initial  successful  demonstration  of 
these  phenomena  in  solid  films  of  DNA,  RNA,  short-chained  oligonucleotides  of  known 
base-pair  sequences  and  spores  was  extended  to  include  proteins,  cells  and  tissues  and  bio 
materials  in  liquid  phase. 

•  Reflection  of  double-  and  single-stranded  DNA  from  herring  and  salmon  as  well  as  samples 
of  protein  Ovalbumin  and  Bacillus  Subtillus  spores  were  measured  and  dielectric  properties 
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of  biological  molecules  in  the  Terahertz  gap  (refractive  index  and  absorption  coefficient 
spectra)  were  characterized. 

•  First  measurements  of  BG  cells  in  the  form  of  aerosol  were  conducted  using  FT 
spectroscopy  which  demonstrated  very  high  sensitivity  for  low  material  concentration. 
Results  of  characterization  were  transfer  to  the  Goodrich  Co  for  use  in  developing  a  system 
for  remote  sensing. 

•  Application  of  TFIz  Spectroscopy  for  biomedical  research:  characterization  of  cancer  cells 
and  tissue. 

•  Characterization  of  Proteins  to  demonstrate  the  capability  of  TFIz  FT  spectroscopy  of 
discriminating  between  different  protein  conformations  and  of  monitoring  the  folding¬ 
unfolding  process. 

•  Developing  TFIz  characterization  technique  for  diluted  solutions  of  biological  materials  for 
possible  future  applications  in  real  time  monitoring  of  biological  processes  on  molecular 
level  and  detecting  water  contaminations. 

•  Further  development  of  theoretical  modeling:  developing  realistic  modeling  approach  and 
application  of  Amber  and  other  available  packages  to  calculate  the  spectra  of  short-  and 
macromolecules. 


•  The  new  TFIz  characterization  technique  for  diluted  solutions  of  biological  materials, 
developed  a  new  technique  for  structural  characterization  of  DNA,  proteins  and  related 
materials  in  diluted  solutions  by  taking  advantages  from  lower  liquid  water  absorption  at  the 
lower  frequency  end  of  TFIz  vs.  IR  and  far  IR  regions.  Reduced  the  amount  of  biomaterial 
in  one  sample  required  for  characterization  or  detection.  The  simple  technique  for  liquid 
sample  preparation  was  developed:  samples  with  less  than  10  ml  of  water  solutions  kept 
between  two  polyethylene  (PE)  or  polycarbonate  (PC)  films  ~10  mm  thickness,  with  the 
concentration  from  0.01  mg/ml  or  0.001%.  The  amount  of  solution  is  controlled  by  the 
spacer  with  the  thickness  from  12  mm. 


Scientific  Personnel  Supported  by  Project: 


Faculty 

Dr.  Tatiana  Globus 
Dr.  Boris  Gelmont 
Dr.  Alexei  Bykhovski 
Dr.  Maria  Bykhovskaia 
Dr.  Tatyana  Khromova 


Graduate  Research  Assistants 
Xiaowei  Li 
Ying  Luo 

Ramakrishnan  Parthasarathy 
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Publication  List : 


•  X.  Li,  A.  Bykhovski,  B.  Gelmont,  T.  Globus,  D.Woolard,  M.  Bykhovskaia,  “Computational  Methods 
for  Analysis  of  Terahertz  Spectral  Signatures  of  Nucleic  Acids”,  5th  IEEE  Conference  on  Nanotechnology, 
Nagoya,  Japan,  July  11-15,  2005,  IEEE  PRESS  BOOKS,  vol.  1, 221  -  224  (2005). 

•  A.  Bykhovski,  X.  Ei ,  T.  Globus,  T.  Khromova  ,  B.  Gelmont ,  D.  Woolard  ,  A.  C.  Samuels,  and  J. 

O.  Jensen,  “THz  absorption  signature  detection  of  genetic  material  of  E.  coli  and  B.  subtilis",  SPIE 
Proceed.  V  5995,  pp.59950N-l:10,  2005,  Optics  East,  Boston,  Oct.  2005. 

•  R.  Parthasarathy,  T.  Globus,  T.  Khromova,  N.  Swami,  D.  Woolard,  “Dielectric  properties  of 
biological  molecules  in  the  Terahertz  gap”,  J.Appl.  Phys.  Lett.,  v.  87,  11,  2005. 

•  T.Globus,D.Woolard,T.W.Crowe,  T.Khromova,  B.  Gelmont,  J.  Hessler,  “Terahertz  Fourier  transform 
characterization  of  biological  materials  in  a  liquid  phase”,  J.  Phys.  D:  Appl.  Phys.  39  3405-3413  (2006). 

•  E.R.  Brown,T.B.  Khromova,  T.  Globus,  D.L.  Woolard,  J.O.  Jensen  and  A.  Majewski,  “THz- 
Regime  Attenuation  Signatures  in  Bacillus  subtilis  and  a  Model  Based  on  Surface  Polariton  Effects”,  IEEE 
Sensors  Journal,  Febr.  2006. 

•  T.  Globus,  T.  Khromova  ,  B.  Gelmont ,  D.  Woolard  ,  and  E.  K.  Tamm,  “Terahertz  characterization  of 
diluted  DNA  solutions”.  Proceedings  of  SPIE  Photonics  West  2006  “Biomedical  Vibrational 
Spectroscopy  III:  Advances  in  Research  and  Industry”,  v.  6093,  Biomedical  Optics  Symposium,  San 
Jose,  California,  21-26  January  2006. 

•  T.Globus,  T.  Khromova,  D.  Woolard  and  A.  Samuels,  “THz  Resonance  Spectra  of  Bacillus 
Subtilis  Cells  and  Spores  in  PE  Pellets  and  Dilute  Water  Solutions,”  SPIE  Proc.  V.62 12-21,  Defense  and 
Security  Symposium  (DSS06);  Conference:  Terahertz  for  Military  and  Security  Applications  IV,  Orlando, 
Florida,  Apr.  17-21,  2006. 

•  T.  Globus,  T.  Khromova  ,  D.  Woolard  and  A.  Samuels,  “THz  Resonance  Spectra  of  Bacillus 
Subtilis  Cells  and  Spores  in  PE  Pellets  and  Dilute  Water  Solutions”,  SPIE  Defense  and  Security 
Symposium,  V  6212-21,  Orlando,  Florida,  17-21  April  2006. 

•  A.  Bykhovski,  T.  Globus,  T.  Khromova  ,  B.  Gelmont ,  D.  Woolard  ,  and  M.  Bykhovskaia  ,  “An 
Analysis  of  the  THz  Frequency  Signatures  in  the  Cellular  Components  of  Biological  Agents”,  SPIE 
Defense  and  Security  Symposium  2006,  V  6212-18,  p.  132-141.  Orlando,  Florida,  April  2006. 

•  D.  E.  Woolard,  Y.  Euo,  B.  E.  Gelmont,  T.  Globus  and  J.  O.  Jensen,  “Bio-molecular  inspired 
electronic  architectures  for  enhanced  sensing  of  thz-frequency  bio-signatures”.  International  Journal  of 
High  Speed  Electronics  and  Systems  (IJHSES),  Vol.20,  No.  10,  1-15  (2005). 

Conclusions  : 

•  Demonstration  of  the  fundamental  character  of  resonance  features  in  transmission 
(absorption)  spectra  of  biological  materials  in  the  THz  spectral  range. 

•  These  resonances  arise  from  electromagnetic  wave  interactions  with  the  low 
frequency  vibrational  modes  within  biological  macromolecules  involving  the 
weakest  hydrogen  bonds  and/or  non-bonded  interactions  between  different 
functional  groups  or  even  between  molecules. 

•  THz  transmission  (absorption)  spectroscopy  of  biological  macromolecules  can 
reveal  these  multiple  resonances  and  the  spectra  can  be  used  for  detection  and  for 
characterization  of  bio-polymers. 
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•  Demonstration  that  Fourier  Transform  spectroscopy  can  be  used  to  measure  low 
frequency  vibrational  modes  of  biomolecules  in  dilute  water  solutions  at  the  THz 
frequencies  10-25  cm  -1.  SE)ur  results  indicate  the  interaction  between  biological 
molecules  and  substrate  material  is  significant  and  likely  controlled  by  substrate 
surface  charges.  This  interaction  causes  the  first  monolayer  of  biomolecules  to 
adsorb  to  the  substrate  in  specific  orientations. 

•  Orientation  effects  contribute  significantly  to  the  resonant  features  in  the  THz 
spectra. 

•  The  band  around  14  cm'*  is  very  sensitive  to  the  orientation  of  the  DNA 
molecules  and  can  be  used  as  a  key  to  control  the  preferable  direction  of  material 
orientation. 


Freqjcrcj.  cm' 
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^smu'ation  result:  usrnf  Jumna  t^reeni  and  Amber 
'bSuei. 
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Good  reproduabiliry  of  spectral  features  of  a 
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Dtfforen:  concentration:  and  a\  erased 
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Solid-State  Transit  Time  THz  Sources 
X.  Bi ,  J.  East  and  G.  Haddad 
Solid  State  Electronics  Eaboratory 
The  University  of  Michigan 

This  report  is  a  brief  summary  of  the  results  of  the  project.  The  details  of  the  research 
are  described  in  a  Ph.D.  Thesis,  “Analysis  and  Design  of  Si/SiGe  Terahertz  Diodes”  by 
Xiaochuan  Bi,  a  student  that  was  supported  by  the  project. 

Summary: 

The  goal  of  the  research  in  this  project  was  to  investigate  novel  THz  sources.  The 
effort  started  with  an  investigation  of  vacuum  Thz  transit  time  diodes.  These  devices 
overcome  many  of  the  limitations  of  semiconductor  based  transit  time  device  oscillators. 
In  particular  the  drift  region  dielectric  constant  is  the  vacuum  value  and  the  drift 
velocities  are  ballistic  rather  than  scattering  limited.  The  lower  dielectric  constant  will 
lower  the  capacitance  per  unit  area  and  the  higher  velocity  associated  with  ballistic 
transport  will  increase  the  drift  region  length,  further  reducing  the  cold  capacitance.  This 
combination  will  allow  a  larger  device  area  for  a  given  impedance  level  and  frequency, 
which  should  increase  the  power  available  from  the  device.  The  ballistic  nature  of  the 
electron  transport  is  also  an  advantage.  In  a  conventional  semiconductor  based  device  the 
carrier  velocity  is  near  the  saturated  value  and  the  resulting  induced  current  is 
approximately  rectangular.  In  the  vacuum  structure  the  acceleration  rather  than  the 
velocity  is  nearly  constant  and  the  resulting  induced  current  is  triangular,  improving  the 
power  generation  in  the  second  14  of  the  RE  cycle.  Detailed  numerical  results  indicated 
that  this  new  device  could  produce  useful  power  at  Thz  frequencies.  However,  we  had 
major  problems  when  we  tried  to  fabricate  experimental  structures.  The  critical  part  of 
the  experimental  structure  was  the  field  emitter  cathode.  This  cathode  was  based  on  a 
micromachined  silicon  structure  that  used  nanoscale  silicon  tips  as  field  emitters.  Earlier 
research  showed  that  these  tip  arrays  could  provide  modest  current  densities.  The  earlier 
work  was  based  on  using  the  emitters  for  very  low  frequency  display  applications. 
However,  although  some  of  the  experimental  structures  produced  current  densities  high 
enough  to  be  useful  the  experimental  results  were  not  reproduceable  and  possible 
improvements  to  the  structure  would  introduce  additional  capacitance  that  would  short 
out  the  RE  portion  of  the  transit  time  structure.  Therefore,  we  decided  to  refocus  the 
research  on  more  conventional  semiconductor  based  structures.  These  results  will  be 
described  next. 

Transit  time  sources  have  a  long  history  but  there  are  several  reasons  to  revisit  them  for 
use  as  very  high  frequency  sources.  Early  devices  were  based  on  silicon  technology  but 
recent  research  has  focused  on  III-V  materials,  partly  because  of  the  use  of  MBE  to 
realize  very  precise  doping  profile  control.  However  improvements  in  silicon  MBE 
associated  with  improvements  in  SiGe  HBT’s  overcomes  this  disadvantage.  The  material 
properties  of  the  III-V  devices  were  better  for  microwave  frequencies.  However,  these 
advantages  can  become  limitations  at  THz  frequencies.  Einally,  at  higher  frequencies  a 
more  detailed  physical  model  is  needed  to  describe  the  transport.  An  electron  traveling  at 
a  velocity  of  lO^cm/sec  in  an  electric  field  of  1  MV/cm  will  gain  1  eV  of  energy  in  10 
nanometers  and  0.1  psec.  These  distance  and  time  scales  can  be  neglected  at  microwave 
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and  millimeter  wave  frequeneies  but  beeome  important  approaehing  THz  frequeneies.  A 
major  portion  of  the  effort  in  the  projeet  involved  a  detailed  transport  level  investigation 
of  the  limitations  of  THz  transit  time  deviees. 

The  researeh  implemented  a  Monte  Carlo  simulation  of  Si  and  SiGe  transit  time  deviees. 
The  results  showed  two  important  transport  effeets  that  are  not  aeeurately  deseribed  by 
simpler  teehniques.  The  first  is  veloeity  overshoot.  Veloeity  overshoot  improves  the 
performanee  of  most  transistors  but  degrades  transit  time  deviee  operation  by  inereasing 
the  indueed  eurrent  near  the  injeetion  phase.  The  design  goal  beeomes  to  reduee  the 
transient  overshoot  as  mueh  as  possible.  This  means  that  silieon  with  a  modest  veloeity 
overshoot  is  a  better  ehoiee  for  THz  transit  time  deviees  than  III-V  materials.  Overshoot 
ean  also  be  redueed  by  injeeting  earners  at  a  higher  energy.  This  means  that  higher 
injeetion  energies  after  ionization  is  better  than  lower  injeetion  energies  after  tunneling 
injeetion.  Both  of  these  results  are  new. 

The  seeond  important  transport  effeet  is  energy  eonservation  or  “dead  zones”.  At  lower 
frequeneies  the  generation  region  ean  be  defined  with  a  high  field  using  the  doping 
profile.  However,  in  higher  frequeney  deviees  the  voltages  are  lower  and  the  ability  to 
define  the  generation  region  with  the  doping  profile  is  redueed  as  the  generation  spills 
into  the  drift  portion  of  the  deviee.  Finally  for  very  short  struetures  the  generation  extends 
over  most  of  the  deviee.  Again  a  Monte  Carlo  analysis  that  ineludes  energy  eonservation 
effeets  is  needed  to  aeeurately  model  deviee  performanee.  A  detailed  analysis  of  silieon 
deviees  operating  between  200  and  500  GHz  was  earned  out  using  the  Monte  Carlo 
simulation.  It  was  clear  that  generation  region  design  limitations  associated  with  energy 
conservation  were  the  fundamental  limit  on  extending  the  high  frequency  operation.  At 
this  point  we  investigated  an  alternative  generation  region  design,  a  SiGe  heterostructure. 
This  design  will  have  a  lower  ionization  rate  material  between  two  higher  ionization  Si 
drift  regions.  This  structure  will  have  a  better  carrier  generation  confinement  and  better 
induced  current  waveform.  A  variety  of  SiGe  based  structures  were  analyzed  to  optimize 
the  operation  above  200  GHz.  The  large-signal  results  showed  that  the  SiGe  MITATT 
diode  improves  the  power  generation  by  63%  at  200  GHz  and  20%  at  250  GHz.  For 
frequencies  beyond  300  GHz,  both  Si  and  SiGe  MITATT  diodes  showed  great  promise 
for  power  generation.  They  produce  more  than  8  mW  at  300  GHz  and  more  than  I  mW  at 
400  GHz. 

List  of  Publications  Resulting  from  This  Work: 

“Analysis  and  Design  of  Si  Terahertz  Transit  Time  Devices”,  X.  Bi,  J.  East  U. 

Ravaioli  and  G.  Haddad,  Solid  State  Electronics,  Issue  5,  May  2006,  pages  889-896. 

“Analysis  and  Design  of  Terahertz  Transit  Time  Diodes,”  X.  Bi,J.  East,  U.  Ravaioli 
and  G.  Haddad,  16*  International  Conference  on  Space  Terahertz  Technology,  May 
2-4,  2005,  Goteborg,  Sweden. 

“Analysis  and  Design  of  Si/SiGe  Terahertz  Diodes”,  Xichuan  Bi,  Ph.D.  thesis.  The 
University  of  Michigan,  2007. 
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Science  and  Technology  of  Chemical  and  Biological  Sensing  at  THz  Frequencies 


Final  Report  Contribution 
Principal  Investigator:  E.R.  Brown,  Ph.D. 

University  of  California,  Santa  Barbara  93106 

Summary 

Three  highlight  areas  were  ehosen  so  summarize  the  researeh  findings  made  by 
the  group  of  E.R.  Brown,  first  at  UCEA  between  2001  and  2004  and  then  at  UCSB 
between  2004  and  2006.  These  are:  (1)  the  development  of  the  photomixing  eoherent 
transeeiver  and  its  applieation  to  THz  speetroseopy,  (2)  high-resolution  detection  and 
modeling  of  spores  of  Bacillus  subtilis,  and  (3)  detection  of  surprisingly  narrow  THz 
absorption  resonances  in  polysaeeharides,  partieularly  lactose  monohydrate. 

1.  THz  Photomixer  Coherent  Transceiver 

A  spectrometer  needs  both  a  souree  to  illuminate  a  sample  material  of  interest  and 
a  deteetor  whieh  measures  how  mueh  gets  through  the  sample.  In  this  projeet  we  used 
two  ErAs:GaAs  photoeonduetive  mixers  for  the  souree  and  deteetor  (Eig.  I).  The  system 
is  driven  by  two  NIR  lasers  whieh  makes  it  highly  eoherent  (coherenee  bandwidth  is  <20 
MHz),  and  therefore  eapable  of  measuring  extremely  narrow  speetral  features.  It  is 
tunable  from  DC  to  3  THz  within  seeonds.  The  two  laser  beams  are  offset  in  wavelength 
and  mixed  together  in  a  biased  photoeonduetors  to  generate  a  cw  THz  beam  at  the  lasers' 
differenee  frequeney.  This  is  similar  to  differenee  frequency  generation  in  nonlinear 
erystals  but  since  bias  is  applied  to  the  photoconduetor,  the  mixing  gains  power  from  the 
applied  bias.  Also,  sinee  the  mixing  oeeurs  within  a  short  length  in  the  semieonduetor  it 
eliminates  the  requirement  for  phase  matehing.  Due  to  the  reeiproeity  theory  in 
eleetromagneties,  the  role  of  the  bias  and  the  THz  beam  in  a  photomixer  in 
interehangeable.  This  can  be  seen  in  Eig  1.,  where  the  transmitting  photomixer,  generates 
THz  power  when  biased,  whereas  the  reeeiver  (Rx)  on  the  right  side,  generates  eurrent 
when  THz  power  is  incident.  Ultrafast,  low-temperature-grown  (ETC)  GaAs 
photoeonduetors  are  now  commonly  used  for  making  short  broadband  pulses  in  switehes 
for  time  domain  speetroseopy  systems  (TDS)  but  in  this  system  we  use  a  superior 
ErAs:GaAs  material  system  by  eollaboration  with  the  material  seienee  department  at 
UCSB.  The  advantage  of  the  transeeiver  is  that  one  ean  get  the  amplitude  as  well  as  the 
phase  response  of  the  material.  To  get  information  about  the  phase  ehange  in  the  system, 
a  retromirror  is  moved  to  ehange  the  optieal  path  on  the  transmitting  photomixer  by  a  few 
THz  wavelengths.  This  renders  a  sinusoidal  interferogram.  With  suffieient  signal  strength, 
the  frequeney,  phase  and  amplitude  of  the  wave  can  be  determined.  However,  sinee  the 
frequeney  ean  be  easily  determined  with  NIR  analysis  of  the  pump  laser  beams,  this  ean 
be  used  to  ereate  a  matehed  filter  that  eorrelates  the  sinusoidal  interferogram  with  the 
interferogram  to  be  expeeted.  This  enables  phase  and  amplitude  determination  even  when 
the  noise  is  larger  than  the  signal. 

Eig.  2  shows  the  spectral  response  if  the  transceiver  during  a  fast  sean  with  broad 
frequency  coverage.  The  maximum  eurrent  in  the  reeeiver  is  5  nA  when  the  THz  power  is 
approximately  10  mierowatts.  The  statie  (iRx/Popt  )  and  differential  sensitivities 
(diRx/dPopt)  are  20pA/W^^^  and  40  mA/W  respectively.  As  a  spectroscopieally  qualitative 
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test  for  the  performance  of  the  transceiver  system,  the  transmission  through  air  between 
1.2  and  1.8  THz  was  measured  as  shown  in  Fig  3.  This  frequency  region  was  chosen 
because  it  contains  well  understood  water  vapour  absorption  lines.  At  each  frequency 
data  point,  the  interferogram  was  measured  while  dithering  a  retro-mirror  by  2  cm  over  a 
10s  period  at  each  frequency  point.  Fig.  3  shows  that  the  amplitude  data  (blue  dots)  agree 
very  well  with  theoretical  model  based  on  HITRAN  parameters  (solid  black  line),  both  in 
frequency  location  and  intensity.  The  ratio  maximum  signal  measured  to  the  minimum  is 
over  50dB.  On  the  right  side  figure  the  change  in  index  of  refraction  of  the  air  in  the  path 
(blue  dots)  is  plotted  along  with  theory  (solid  red  line).  This  shows  that  the  system  is 
successfully  acquiring  the  phase,  tracking  even  rapid  changes,  almost  discontinuities  in 
the  index  of  refraction. 
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FIG.  1.  Experimental  setup.  The  TFIz  beam  is  generated  by  transmitter  (Tx),  penetrates 
sample  and  is  converted  to  current  on  the  receiver  (Rx). 


FIG.  2.  Spectral  response  of  coherent  transceiver. 
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Fig  3.  Amplitude  (left)  and  phase  (right)  response  of  transeeiver  when  seanning  aeross 
atmospherie  water  vapor  lines.  The  solid  line  is  a  Lorentzian  model  based  on  parameters 

from  the  HITRAN  database. 

II.  Measurement  and  Modeling  of  Narrow  Resonances  in  Bacillus  Subtilis 

Mueh  of  our  research  during  the  MURI  concentrated  on  bioparticles,  especially 
bacterial  spores  of  the  genus  Bacillus.  Scanning  electron  microphotographs  (SEMs)  of  a 
typical  spore  of  B.  subtilis  var.  niger  (formerly  known  as  Bacillus  globigii),  or  BG  for 
short,  indicate  that  it  has  a  spheroidal  shape  with  a  length  of  «  1.2  pm  and  width  of  «  0.8 
pm.’  The  spore  outer  coat  (OC)  represents  approximately  30  to  60%  of  the  dry  weight  of 
the  spore.  The  OC  is  a  tough  proteinaceous  layer  having  an  unusually  high  content  of 
cysteine  and  other  hydrophobic  amino  acids.  This  makes  the  OC  highly  resistant  to 
protein  solvents,  and  also  contributes  to  its  rather  high  atomic  density  and  dark 
appearance  in  SEMs.  Also  important  to  the  present  experiments  is  the  reported 
crystallinity  of  this  layer.  As  in  other  solids,  crystallinity  tends  to  support  collective 
excitations,  such  as  the  phonons,  and  renders  them  with  much  longer  lifetime  than  occurs 
in  polycrystalline  or  ceramic  materials. 

In  the  MURI  research  samples  of  BG  were  prepared  using  spores  mixed  with 
vegetative-cell  material  and  packaged  in  three  different  forms:  (1)  BG  material  only 
(vegetative  cell  matter  plus  spores),  (2)  BG  material  distributed  in  polyethylene  (PE) 
powder,  and  (3)  BG  material  suspended  in  air  as  an  aerosol  (done  at  the  University  of 
Virginia,  Dr.  T.  Globus).  The  most  important  data  was  obtained  on  the  diluted  BG 
sample.  Estimates  of  the  mass  fraction  of  this  sample  and  the  two  others  are  listed  in 
Table  I.  Optical  microscopic  inspection  of  the  BG  material  showed  that  the  vegetative¬ 
cell  material  were  primarily  in  the  form  of  pieces  of  broken  cell-wall  remnants. 

For  the  purpose  of  extracting  attenuation  signatures,  the  desired  quantity  from  the 
measurements  was  the  transmission  spectrum  through  each  BG  sample  normalized  to  the 
transmission  spectrum  through  its  control  sample.  In  principle,  such  normalization 
eliminates  the  effects  of  frequency-dependent  source  power,  sample-dependent  reflection 
coefficient,  etc.  The  normalized  transmission  T’(v)  through  the  dilute-BG  sample  is 
plotted  in  Fig.  4  between  200  GHz  and  just  under  1.08  THz,  above  which  water  vapor 
absorption  reduced  the  SNR  to  unacceptable  levels.  The  frequency  registration  and 
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normalization  procedure  does  a  good  job  of  eaneeling  the  water  absorption,  but  small 
undulations  in  T’(v)  still  oeeur  aeross  the  speetrum.  These  are  attributable  at  least  in  part 
to  small  ehanges  in  the  instrumental  standing  wave  pattern  eaused  by  the  ehanging  of 
samples  and  the  slight  difference  in  their  plaeement  and  refractive  index.  Provided  that 
the  refraetive  index  and  position  is  not  ehanging  with  frequency,  this  type  of  undulation 
is  expeeted  to  lie  inside  of  an  envelope  defined  roughly  by  the  two  dashed  eurves  in  Fig. 

4.  The  vertieal  range  of  the  envelope  is  «0.8  dB  (17%) 

Three  attenuation  signatures  fall  -'10%  below  the  lower  extent  of  the  envelope  as 
highlighted  by  the  dotted  reetangles.  They  are  eentered  around  250,  415,  and  1035  GHz, 
respeetively,  are  progressively  strong  in  attenuation,  and  have  all  been  reported  by  other 
researchers,  particularly  Dr.  T.  Globus  of  U.  Virginia.  The  first  two  signatures  are 
regarded  as  probable  deteetions  sinee  they  eorrespond  to  the  275  and  425  GHz  signatures 
of  the  eoneentrated-BG  sample.  The  1.035  THz  signature  is  regarded  as  the  most 
signifieant  sinee  it  is  stronger  than  the  other  two  and  has  also  been  measured  recently  by 
time-domain  THz  spectroscopy."  Again,  AT’  is  derived  by  dividing  out  the  broadband 
attenuation,  this  time  taken  as  an  average  through  to  the  eenter  of  the  envelope.  The 
eorresponding  spore  eoneentration  and  attenuation  eross  seetion  are  listed  in  Table  I.  In 
eomparing  the  dilute  and  eoneentrated-BG  results,  it  is  remarkable  that  10  times  less  BG 
material  (by  mass)  is  produeing  3  times  stronger  attenuation  per  spore.  The  linewidths 
are  diffieult  to  determine  preeisely,  so  are  merely  bounded  in  Table  1  as  <  15  GHz  (250 
GHz  signature)  and  <10  GHz  (415  GHz  signature).  They  are  somewhat  narrower  than 
the  corresponding  linewidths  in  the  concentrated  sample. 

Additional  work  done  under  the  MURI  developed  a  physieal  model  for  the  origin 
of  the  signatures  observed  in  BG  spores.  The  model  is  based  on  the  existenee  of  polar 
transverse  optieal  (TO)  phonons  in  the  polyerystalline  outer  eoat  (OC)  of  eaeh  BG  spore. 

The  interaetion  of  these  TO  phonons  with  THz  radiation  was  analyzed  by  invoking 
conservation  of  energy  and  momentum,  leading  to  a  strong  eoupling  effeet  between  them 
ealled  a polariton.  The  resulting  dispersion  relation,  cos  vs  K  (K  being  the  polariton 
wave  veetor  in  the  biopartiele),  is  typieal  of  eoupled-mode  systems  having  two  branehes, 
and  is  shown  for  our  BG  spore  model  in  Fig.  5(a).  Strong  eoupling  oecurs  only  in  and 
around  the  “gap”  between  the  branches,  the  lower  extent  being  cot,  and  the  upper  extent 
being  «[(cot)  +  (cOp)  ]  .  An  important  aspeet  of  the  polaritonie  eoupling  is  the  small 
size  of  the  gap  eompared  to  the  cot.  This  is  eaused  by  the  small  value  of  cOp  or, 
equivalently,  the  large  mass  of  the  unit  eell.  This  has  one  important  eonsequenee  on  the 
THz  op  ties:  the  range  of possible  polariton  frequencies  associated  with  a  particular  TO- 
phonon  resonance  is  very  small. 

The  overall  effect  on  the  propagating  THz  radiation  oeeurs  through  the  dieleetrie 
funetion  whieh  for  polaritons  and  other  harmonie-oseillator-lrke  exeitations  in  solids  is 
often  approximated  in  a  first  analysis  by  the  Lorentzian  expression 

s(co)  =  So  +  — - ^ ^ — 

0)^-0}  -  jyo) 

where  8o  is  the  baekground  dieleetrie  eonstant  in  the  limit  of  high  frequeneies  (»  cot),  y  is  the 

2  1/2 

damping  rate,  cop  =(pcq  /M)  is  the  effeetive  plasma  frequeney,  q  is  the  effeetive  charge  (and 
thus  the  degree  of  polarity),  M  is  the  mass  per  primitive  unit  eell,  and  pc  is  the  density  of  unit 
eells  and  hence  of  TO  phonons.  In  bulk  solids  this  would  suffiee  in  calculating  the  THz 
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attenuation  properties.  But  in  the  present  samples  we  must  also  aecount  for  the  particulate 
nature  of  the  BG  spores,  whieh  is  greatly  faeilitated  by  three  facts.  The  first  is  that  the  spore  size 
is  ~300  times  less  than  the  radiative  wavelengths  of  interest  in  the  sub-THz  region.  This  allows 
us  to  evaluate  the  electromagnetic  interaction  in  the  “electrostatic”  limit  whereby  the  exciting 
and  scattered  electric  fields  (E)  are  assumed  to  be  spatially  uniform  around  the  particle.  The 
second  fact  is  that  the  spore  shape  is  spheroidal,  which  from  electrostatics  makes  the  polarization 
field  induced  in  the  spore  uniform!"  The  third  fact  is  that  the  spores  in  our  dilute  and  aerosol 
samples  are  rarefied  enough  that  they  can  be  treated  as  independent  scatterers.  The  spores  then 
contribute  to  a  radiative  extinction  coefficient.  A,  simply  by  A  =  PvCa  where  pyis  their 
volumetric  density  and  Ca  is  the  absorption  cross  section  plotted  vs  frequency  in  Fig.  5(b). 
Further  details  of  the  model  can  be  found  in  the  paper:  “THz  Regime  Attenuation  Signatures  in 
Bacillus  subtilis  and  a  Model  Based  on  Surface  Polariton  Effects,”  E.R.  Brown  et  al,  IEEE  Sensors  J., 
Vol.  6,  No.  5,  pp.  1076-1083  (2006). 


Table  I. 


Spectral 

Feature 

[GHz] 

BG 

Sample 

Estimated 

Spore 

Density  [cm'^ 

AT’ 

Ca” 

[cm^] 

Av 

[GHz] 

concentrated 

-7x10“ 

0.40 

2 . 2x10'"" 

20 

-260 

dilute 

-7xl0' 

0 . 75 

6.8x10'"" 

<  15 

concentrated 

-7x10“ 

0.30 

2 . 9x10'"" 

»  10 

-420 

dilute 

-7xl0' 

0.70 

8.5x10'"" 

<  10 

*  measured  at  minimum  transmission  of  resonant  signature 
**estimated  attenuation  cross  section  per  spore. 


Fin.  4. 
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Fig.  5. 

Polysaccharides 

Most  of  the  biomass  on  the  planet  consists  of  saccharides.  Understanding  their 
THz  dynamics  should  therefore  be  of  high  significance  to  bimolecular  detection.  IR 
spectroscopy  is  used  to  examine  the  vibrational  absorption  spectrum  of  organic  molecules. 
The  frequency  of  these  transitions  is  related  to  the  strength  of  the  chemical  bond  and 
mass  of  the  bonded  atoms.  At  lower  frequencies,  and  specifically  in  the  THz  regime,  one 
expects  rotational  spectral  features  or  vibrational  spectra  of  weaker  bonds  and  larger 
masses.  Polysaccharides  have  attracted  attention  in  THz  spectroscopy  because  they 
consist  of  hydrogen  bonded  networks  of  sugar  molecules,  saccharides.  Each  sugar 
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molecule  is  composed  of  two  or  more  monosaccharides  bonded  together  with  an  oxygen 
molecule.  This  gives  rise  to  a  large  degree  of  freedom  for  oscillation  within  the  sugar 
molecule. 

One  of  the  interesting  things  about  sugars  is  that  a  very  large  family  of  anomers 
has  exactly  the  same  chemical  composition  but  the  alignment  of  chemical  bonds  makes 
the  molecular  dynamics  and  the  biological  role  completely  different.  Because  of  the  scale 
of  the  molecules  THz  spectroscopy  is  contains  the  ideal  frequency  range  to  directly  detect 
the  molecular  dynamics  and  by  comparing  the  spectra  of  different  sugars. 

We  measured  an  absorption  line  in  lactose  around  518  GHz  with  a  photomixing  source 
and  a  room  temperature  bolometer  as  a  detector.  The  samples  were  prepared  by  mixing 
10%  lactose  and  90%  polyethylene  powder  (PE)  in  a  10mm  thick  pellet.  The  spectrum  is 
shown  in  Fig.  6.  The  line  shape  is  fitted  with  a  Lorentzian  peak: 


Liu)  = 


J  D 

;r  +  {o-Of^Y 


And  the  fitting  parameters  for  the  half-width,  D,  were  16  GHz,  the  total  intensity,  A  was 
0.1 12  and  the  center  was  at  533  GHz.  Fig.  6  also  shows  a  spectrum  taken  without  sample 
in  the  path  and  the  residual  error  in  the  Forentzian  fitting.  The  linewidth  of  sub-THz 
absorption  features  (generally  from  lattice  vibrations)  in  solids  is  excepted  to  inversely 
proportional  to  the  scattering  time  of  optical  phonons.  The  line  width  of  absorption 
features  might  thus  have  implications  on  the  solid-state  physics  of  the  material,  in 
particular,  the  interaction  of  phonons  and  the  phonon  density  of  states.  In  this  case,  the 
life  time  of  this  transition  is  1/(27t:  -IhGHz)  ~  70  ps. 

Fig.  7  shows  the  extraction  of  three  absorption  lines  in  Maltose  monohydrate 
using  a  Forentzian  fit  and  assuming  a  Mie  scattering  background  trend.  Mie  scattering  is 
often  seen  in  infrared  spectroscopy  of  aerosol  sub-wavelength  particles.  In  our  case  we 
have  higher  dielectric  powdered  saccharides  in  a  lower  dielectric  polyethylene  powder 
matrix  which  has  similarly  Mie  behaviour.  The  three  Forentzian  peaks  that  were 
numerically  extracted  in  maltose  were  located  at:  1074,  1593,  and  1985  GHz. 

The  contrast  between  these  two  samples,  lactose  and  maltose  which  have  the  same 
composition,  but  completely  different  THz  spectra,  demonstrate  how  rich  the  THz  spectra 
are  in  information  about  molecular  structure  beyond  the  composition. 
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Fig,  6  Absorption  in  a-lactose  monohydrate.  The  main  graph  shows  the  power  spectrum 
of  a  photomixer  illuminating  a  cryogenic  bolometer  with  a  lactose  sample  in  between. 
The  two  insets  show  a  fit  to  theoretical  Lorentzian  functions  at  0.533THz  and  1.378  THz. 
There  is  also  possibly  an  absorption  feature  right  below  1 .2THz. _ 
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Fig,  7  Absorption  in  maltose  monohydrate.  The  absorption  spectrum  of  powdered 
maltose  can  be  separated  into  a  superposition  of  three  Lorentzian  absorption  lines  (1074, 
1593,  and  1985  GHz)  and  a  Mie  scattering  background  trend  _ 
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Six-Port  Reflectometers  and  Calibration  Techniques  for  Submillimeter 
Measurement  of  Biological  Materials 


Summary 

The  test  and  measurement  infrastructure  that  has  been  indispensable  to  the  development  and 
characterization  of  RF,  microwave,  and  optical  systems  does  not  exist  for  the  submillimeter-wave 
region  of  the  spectrum.  This  has  long  been  recognized  as  a  major  impediment  to  both  the  use  and 
development  of  this  spectral  region  and  has  impeded  progress  on  the  measurement  and 
characterization  of  materials,  including  chemical  and  biological  agents.  This  research  effort  is 
focused  on  the  development  of  instrumentation,  specifically,  six-port  reflectometers,  to  allow 
direct  measurement  of  complex  scattering  parameters  (magnitude  and  phase)  of  samples  under 
test. 

Waveguide-Based  Reflectometers 

The  most  common  transmission  media  in  the  submillimeter  region  are  rectangular  waveguide  and 
free-space.  Consequently,  the  reflectometers  investigated  under  support  of  this  MURI  were 
designed  using  rectangular  waveguide  as  the  medium.  The  primary  focus  was  on  frequency  bands 
extending  beyond  those  currently  covered  by  commercial  equipment  (300  GHz  and  higher).  A 
photograph  and  detail  of  the  prototype  reflectometer  developed  under  this  grant  is  shown  in  Fig. 
8.  The  instrument  consists  of  a  5  cm  section  of  rectangular  waveguide  and  a  set  of  four  Schottky 
detector  diodes.  Because  only  three  detectors  are  required  to  measure  passive  loads,  three  diodes 
serve  as  the  primary  detectors  and  the  fourth  is  used  as  a  backup.  The  diodes  used  in  this  work  are 
SC1T5  planar  Schottky  diodes  that  were  fabricated  in  the  Microfabrication  Laboratory  at  the 
University  of  Virginia.  These  devices  are  typically  used  for  submillimeter-wave  mixers  operating 
near  600  GHz  and  have  series  resistance  in  the  range  of  6 — 9  Q.  Details  of  the  waveguide  and 
probe  design  are  provided  below: 

Waveguide  Design 

The  waveguide  block  used  for  this  work  was  fabricated  by  Custom  Microwave,  Inc.  The  primary 
waveguide  channel  has  dimensions  of  406  pm  X  787  pm,  resulting  in  single  mode  propagation 
over  the  190  GHz  to  370  GHz  band.  Shallow  (125  pm  deep  X  178  pm  wide)  cross  channels  are 
machined  into  the  block  to  accommodate  the  microstrip  probes  that  sample  the  standing  wave  in 
the  guide.  The  dimensions  of  the  cross  channels  were  chosen  to  suppress  excitation  of  high-order 
modes,  leaving  only  the  quasi-TEM  shielded  microstrip  mode  to  travel  to  the  diode  detectors.  The 
probe  spacing  is  195  pm,  corresponding  to  A,g/6  at  320  GHz.  Moreover,  two  narrow  (100  pm 
wide)  side  channels  are  machined  parallel  to  the  primary  waveguide  channel  (see  Fig.  8(b)). 
These  side  channels  allow  quarter-wave  long  bond  wires  to  be  included  as  dc  returns-to-ground 
for  the  Schottky  power  detectors.  The  waveguide  flanges  are  designed  to  mate  with  standard  WR- 
3  components. 

Probe  Design  and  Assembly 

Sampling  probes  and  microstrip  lowpass  fdters  for  the  reflectometer  were  designed  with  the  use 
of  Ansoffs  High  Frequency  Structure  Simulator  (HFSS).  To  minimize  perturbations  to  the 
standing  wave  in  the  guide  and  loading  by  the  detectors,  the  probes  were  originally  designed  to 
couple  a  maximum  of  5%  of  the  incident  power  in  the  guide  to  the  detectors.  A  second  set  of 
probes  were  also  designed  to  increase  the  coupling  to  15%.  The  probes  and  fdters  were  fabricated 
using  standard  photolithographic  techniques  on  125  pm  thick  quartz  substrates.  Low -pass  fdter 
sections  integrated  with  the  probes  were  designed  to  block  the  270 — 320  GHz  submillimeter 
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signal  while  permitting  dc  biasing  of  the  detectors.  Corrections  for  step  discontinuities  and  the 
effect  of  the  waveguide  cover  shielding  the  microstrip  structure  were  analyzed  and  compensated 
for  with  the  aid  of  HFSS.  The  low-pass  fdter  consists  of  an  10-section  stepped  impedance  line 
with  low  impedance  sections  of  65Q  and  high  impedance  sections  of  145Q. 
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Fig.  8.  (a)  Photograph  of  the  suhmillimeter  $ampled-line  reflectometer.  (h)  Detail  showing  the 
waveguide  dimensions  and  placement  of  the  sampling  prohes. 


(a)  (b) 

Fig.  9.  (a)  Photograph  of  a  completed  detector  circuit  showing  the  flip-chip  mounted  diode, 
waveguide  prohe,  low-pass  filter,  and  dc  return-to-ground,  (h)  Photograph  of  the  reflectometer 
showing  the  positions  of  all  four  sampling  prohes. 


100  pm  wide  gaps  included  in  the  probe  design  accommodate  the  detectors.  Small  amounts  of 
silver  epoxy  are  applied  on  either  side  of  these  gaps  and  the  diodes  are  flip-chip  mounted  across 
the  gaps.  The  epoxy  is  then  cured  in  a  temperature  controlled  oven  at  100  C  for  approximately 
two  hours.  The  final  step  in  assembly  consists  of  mounting  the  probes  to  their  respective  locations 
in  the  block.  This  is  done  using  commercially  available  UV-cured  epoxy.  A  photograph  of  a 
completed  probe  circuit  is  shown  in  Fig  9(a)  and  the  completed  reflectometer  is  shown  in  fig  9(b) 
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Measurements  and  Characterization 


The  experimental  setup  used  to  evaluate  the  sampled-line  re°ectometer  is  shown  in  fig.  10.  A 
commercial  frequency  synthesizer  (model  HP-8350B,  50  MHz — 22  GHz)  is  used  as  the  primary 
source.  This  is  followed  by  a  broadband  frequency  tripler  (model  295-320FEM)  designed  and 
fabricated  by  Virginia  Diodes,  Inc.  An  adjustable  attenuator  placed  before  the  tripler  permits 
control  of  the  input  power  provided  to  the  reflectometer.  DC  bias  is  supplied  to  the  Schottky 
detectors  with  commercial  bias  supplies  (model  E3610A  by  Hewlett  Packard)  and  the  output 
current  of  the  detectors  is  monitored  with  a  set  of  Keithley-2000  6  -  1/2  digit  multimeters 

To  establish  that  the  sliding  termination  used  to  calibrate  the  sampled-line  reflectometer  behaves 
as  desired,  the  power  detected  by  the  Schottky  diodes  was  measured  as  a  function  of  backshort 
position.  Figure  11(a)  shows  the  measured  current  through  one  of  the  detectors  as  a  function  of 
backshort  position.  For  this  measurement  the  input  submillimeter-signal  was  fixed  at  300  GHz.  In 
addition  to  allowing  the  performance  of  the  backshort  to  be  assessed,  the  measurement  shown  in 
fig.  1 1  (a)  shows  there  is  no  significant  harmonic  content  in  the  source  signal  other  than  the 
fundamental.  Because  the  six-port  method  relies  on  direct  detection,  it  can  be  susceptible  to 
spurious  responses  arising  from  the  presence  of  harmonics.  Figure  11(b)  shows  the  measured 
locus  the  sliding  short  traces  out  in  the  reflection  coefficient  (E)  plane  at  300  GHz.  This 
measurement  was  done  using  the  sampled-line  reflectometer  after  it  had  been  fully  calibrated.  For 
comparison,  the  locus  produced  by  the  sliding  load  was  also  measured  using  an  HP8510C  vector 
network  analyzer  outfitted  with  the  Oleson  Microwave  Eaboratory's  (OME)  V03VNA-T/RX2 
module  (shown  in  figure  1 1(b)). 
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Fig.  10.  Experimental  setup  for  testing  and  evaluating  the  sampled-line  reflectometer. 

The  second  calibration  step  for  the  reflectometer  consists  of  using  three  well-characterized 
standards  to  determine  the  complex  coefficients  for  the  standard  one-port  error  model.  In  this 
work,  the  custom-made  short  circuit  and  two  offset  shorts  previously  described  were  used  for  this 
calibration  step. 

Initial  measurements  performed  on  the  system  utlized  the  probes  with  low  coupling  (5%).  Above 
300  GHz,  two  of  the  diode  detectors  exhibited  very  low  responsivities,  near  the  noise  floor  of  our 
measurement  system.  It  is  believed  that  this  behavior  is  likely  a  result  of  misalignment  and 
inaccurate  placement  of  the  diode  chips/probes  associated  with  the  manual  assembly  of  the 
reflectometer.  As  a  result,  scattering  parameter  measurements  in  the  270  GHz  to  300  GHz  range 
are  presented.  Verification  standards  are  not  available  for  the  submillimeter-wave  range. 
Consequently,  operation  of  the  reflectometer  was  evaluated  using  the  offset  shorts  designed  for 
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calibration.  Because  only  two  offset  shorts  (of  lengths  33  mils  and  37  mils)  are  used  in  the 
calibration  procedure  for  the  270  GHz  to  300  GHz  range,  a  third  offset  (70  mils)  can  be  used  as  a 
check  to  verify  the  reflectometer  is  functioning  properly.  Remeasurement  of  the  standards  used 
for  calibration  provides  additional  information  regarding  the  repeatability  of  measurements. 
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Fig.  11.  (a)  Plot  of  the  detected  current  vs.  position  for  the  sliding  load  used  to  calibrate  the 
reflectometer.  For  this  measurement,  the  suhmillimeter  source  was  set  to  270  GHz  and  a  single  diode 
detector  was  used,  (h)  Loci  traced  out  hy  the  sliding  short  in  the  F-plane  measured  with  the  sampled- 
line  reflectometer  and  the  commercially  availahle  OML  V03VNA-T/RX2  module 
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Fig.  12  (a)  Measured  phase  vs.  frequency  for  three  offset  short  circuits,  (h)  Measured  magnitude  vs. 
frequency  for  three  offset  short  circuits.  Also  shown  is  the  0  dB  reference  for  lossless  terminations. 

Figure  12(a)  shows  the  measured  phase  of  si  1  for  the  three  offset  shorts  over  the  270  GHz — 285 
GHz  band.  The  measured  phase  is  close  to  that  predicted  using  an  ideal,  lossless  short-circuited 
waveguide  as  a  model  and  follows  the  expected  slope  as  a  function  of  frequency.  It  should  be 
noted  that  no  waveguide  losses  are  accounted  for  in  either  the  theoretical  curves  or  circuit  models 
used  for  the  calibration  standards.  Moreover,  for  these  original  measurements,  a  mechanically- 
tuned  Gunn  diode  source  was  used  to  drive  the  multiplier  chain.  There  is  approximately  a  ±5° 
uncertainty  in  the  phase  measurements  associated  with  use  of  this  Gunn  diode  source  and  this 
uncertainty  arises  from  frequency  drift  of  the  device  as  well  as  manual  imprecision  in  adjusting 
the  micrometers  required  for  frequency  tuning.  Because  the  offset  shorts  used  to  calibrate  the 
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reflectometer  are  ideally  lossless,  the  magnitude  of  the  measured  reflection  coefficient  is  expected 
to  be  near  0  dB.  Magnitude  measurements  performed  on  the  offset  shorts  are  shown  in  fig.  12(b). 
In  general,  the  measurements  for  each  shim  average  near  the  expected  0  dB,  but  some  data 
deviate  substantially,  indicating  significant  measurement  error.  Several  potential  sources  of  error 
exist  in  this  prototype  instrument  and  are  likely  responsible  for  this.  These  include  imprecise 
modeling  of  the  calibration  standards  (which  neglects  loss),  repeatability  in  connecting  the 
waveguide  flanges,  and  the  low  power  levels  detected  by  the  diodes. 

A  Reflectometer  Calibration  Method  Resistant  To  Waveguide  Flange  Misalignment 

To  address  the  calibration  issues  note  above,  a  new  reflectometer  calibration  method  is  described 
that  utilizes  four  standards:  a  flush  short,  two  delay  shorts  with  unspecified  but  different  phases, 
and  an  open-ended  waveguide.  The  calibration  method  eliminates  the  requirement  for  a  precision 
waveguide  matched  load  which  can  be  problematic  to  realize  at  sub-millimeter  wavelengths. 
Importantly,  it  is  shown  that  the  technique  is  resistant  to  waveguide  flange  misalignment,  which 
is  among  the  most  serious  factors  that  degrade  the  calibration  accuracy  of  vector  network 
analyzers  operating  above  100  GHz.  Scaled  measurements  using  this  method  have  been 
performed  in  the  W-band  (75 — 110  GHz)  where  it  is  readily  compared  to  other  calibration 
methods  such  as  thru/reflect/line  (TRL)  to  assess  its  utility  and  performance.  Measurement  results 
demonstrate  the  robustness  of  this  new  calibration  method  and  have  verified  it  is  superior  to  the 
commonly  utilized  short/delay  ed-short/load  (SDL)  technique  with  respect  to  the  influence  of 
waveguide  flange  misalignment.  This  work  was  published  in  the  IEEE  Transactions  on 
Microwave  Theory  and  Techniques  in  July,  2006. 

Vector  network  analyzers  (VNA's)  operating  up  to  the  WR-3  waveguide  band  (220 — 325  GHz) 
are  commercially  available  today.  However,  for  frequency  bands  above  100  GHz  (where 
waveguide  is  the  most  common  transmission  medium)  their  measurement  performance  is  often 
limited  by  the  fabrication  tolerances  of  calibration  standards  that  are  used  along  with  the  precision 
with  which  the  waveguide  flanges  can  be  mated.  Flange  misalignment  (illustrated  in  Fig.  13) 
resulting  from  the  limited  precision  of  waveguide  machining  techniques  is  among  the  most 
serious  challenges  to  repeatable  and  accurate  scattering  parameter  measurements  at  the  upper  end 
of  the  millimeter-wave  spectrum.  In  fact,  waveguide  flange  misalignment  can  often  be  as  large  as 
100  pm  with  standard  machining  and  65  pm  using  precision  fabrication  methods. 


Fig.  13  Illustration  of  waveguide  flange  misalignment  due  to  limited  machining  tolerances  and 
fabrication  errors  associated  with  alignment  pins  and  holes. 
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Flange  misalignment  ean  limit  the  ealibration  aeeuraey  of  VNA's  in  several  ways. 
Foremost,  misalignment  ean  result  in  a  signifieantly  degraded  return  loss  for  waveguide 
matehed  terminations.  As  an  illustrative  example,  HFSS  simulation  has  shown  a  return 
loss  as  low  as  17  dB  of  a  WR-3  waveguide  matehed  load  that  is  misaligned  with  the  test 
port  by  76  pm  along  the  waveguide  E-plane.  Misalignment  of  waveguide  flanges  ean  also 
result  in  eonsiderable  phase  error  in  the  response  of  high-refleotion  loads  (sueh  as  delay 
shorts)  that  are  used  as  ealibration  standards.  HFSS  simulation  has  also  been  performed 
on  a  WR-3  waveguide  quarter-wavelength  offset  short  that  is  misaligned  to  the  test  port 
by  76  pm  along  the  E-plane.  In  this  case  the  phase  error  exceeds  20°.  Einally,  flange 
misalignment  can  also  cause  the  “thru”  and  “line”  standards  used  in  the  “thru/reflect/line” 
(TRE)  calibration  technique  to  be  mismatched  to  the  test  ports,  resulting  in  substantial 
degradation  of  calibration  accuracy. 

Due  to  the  scarcity  of  verification  standards  for  frequencies  above  100  GHz,  the  new 
calibration  technique  was  demonstrated  and  assessed  in  the  W-band  (75 — 110  GHz) 
where  alternative  calibration  methods  such  as  TRE  can  be  used  as  a  basis  for  comparison. 
The  distance  between  the  measured  and  TRL-referenced  data  on  the  E-plane,  AF,  is 
calculated  and  plotted  in  Fig.  14  for  different  calibration  methods.  AF  is  defined  as:  Ar=| 
F-  FtrlI,  where  F  is  the  measured  reflection  coefficient  and  Ftrl  is  the  reflection 
coefficient  corrected  using  the  TRE  calibration  data.  The  robustness  of  the  new  method 
and  its  superiority  to  the  commonly  utilized  SDL  under  the  influence  of  waveguide 
flange  misalignment  have  been  demonstrated. 
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Fig.  14.  AF  (dB)  versus  frequency  for  an  E-H  waveguide  tuner  using  different  calibration  methods. 
The  new  method.  Misaligned  SDL  and  SDDL  are  under  the  influence  of  waveguide  flange 
misalignment.  (The  SDDL  method  was  proposed  by  other  researchers  and  utilizes  a  short,  two  delay 
shorts  and  a  matched  load.) 
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Measurement  Fixtures  for  Calibrated  Measurements  of  Biological  Materials 
Summary: 

Knowledge  of  material  parameters,  such  as  complex  permittivity,  is  often  desired  in  applications 
such  as  chemical  and  biological  spectroscopy.  A  relatively  simple  measurement  scheme  for 
permittivity  was  proposed  by  Nielsen,  who  solved  the  scattering  problem  for  a  dielectric  post 
placed  in  the  center  of  a  waveguide  and  aligned  along  the  E-plane.  Moreover,  an  accurate  and 
general  theory  for  such  a  measurement  scheme  has  been  developed  by  Leviatan.  In  this  work, 
these  methods  are  adapted  to  a  one-port  measurement  system  for  characterizing  biological 
materials  that  makes  use  of  a  calibration  method  designed  to  improve  measurement  accuracy  and 
significantly  reduce  system  complexity. 

The  measurement  scheme  we  adopt  from  Nielsen  is  shown  in  Fig.  15,  where  the  sample 
biological  material  is  placed  inside  a  quartz  tube  with  inner  diameter  of  0.8  mm.  A  waveguide 
short  circuit  is  placed  at  quarter  wavelength  from  the  tube  and  the  reflection  coefficient  of  the 
system  is  measured.  With  proper  calibration,  the  effects  of  the  waveguide  fixture  can  be  removed 
and  the  complex  permittivity  of  the  sample  material  derived  from  the  measurement.  The  test 
fixture  was  fabricated  for  one -port  measurements  in  a  IT-band  waveguide  (75-llOGHz).  The 
constitutive  parameters  of  the  material  in  the  tube  are  determined  by  fitting  the  complex 
permittivity  to  obtain  a  best  fit  between  the  calculated  and  the  measured  reflection  coefficients. 


Typical  one-port  calibration  methods  (such  as  the  short/open/load  technique)  cannot  remove 
errors  associated  with  unrepeatable  or  imperfect  connections  between  the  waveguide  test  fixture 
and  the  network  analyzer  test  port.  Thus,  an  in-situ  calibration  method  is  introduced  in  which  the 
test  fixture  remains  connected  to  the  network  analyzer  test  port  for  both  the  calibration  and 
material  measurement.  As  with  standard  one -port  calibrations,  three  standards  are  chosen;  in  this 
case,  an  empty  waveguide,  an  empty  quartz  tube  inserted  in  the  waveguide,  and  a  metal  (gold- 
coated)  rod  inserted  in  the  waveguide. 
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Fig.  15.  (a)  Photograph  of  the  waveguide  test  fixture  with  quartz  capillary  tube,  (h)  Permittivity  of 
Escherichia  coli  as  fouud  with  the  uew  calihratiou  method. 
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Measurements 


The  calculated  complex  permittivity  of  Escherichia  coli  cell  with  solution  concentration  of  1 
mg/ml  based  on  the  measurement  with  the  new  calibration  procedure  at  25  C  is  given  in  figure 
15(b).  The  validity  of  the  new  calibration  method  was  assessed  by  measuring  deionized  water  as  a 
control.  Comparison  of  the  calibrated  measurements  with  other  available  reference  values  shows 
that  the  measured  permittivities  are  close  with  slight  discrepancies  between  the  imaginary  parts. 
This  is  attributed  to  limitations  in  the  present  fixture  model  which  does  not  account  for  fringing 
fields  that  may  extend  in  the  tube  beyond  the  boundary  of  the  rectangular  waveguide. 


Integrated  Diode  Phase-Shifters  for  Sideband  Generation  at  Suhmillimeter 

Wavelengths 


Background  and  Objective: 

This  effort  was  aimed  at  developing  broadband  tunable  submillimeter-wave  sources  based  on 
sideband  generation.  The  most  efficient  sideband  generators  are  constructed  from  phase-shifters 
capable  of  producing  a  180  degree  phase-shift  of  the  incident  signal.  The  resulting  sidebands 
from  the  phase  modulation  can  be  tuned  over  a  broad  frequency  range  and  used  as  a  scanned 
source  for  spectroscopic  measurements. 

Summary  and  Findings: 

The  phase  shifter  investigated  in  this  study  is  based  on  a  well-known  balanced  architecture 
consisting  of  a  3-dB  90  hybrid  and  a  pair  of  identical  reflective  terminations.  These  terminations 
are  realized  with  integrated  Schottky  barrier  varactors  with  device  parameters  designed  to  yield  a 
reasonable  tradeoff  between  bandwidth,  insertion  loss,  and  amplitude  balance.  The  operating 
band  of  this  prototype  is  chosen  as  170-260  GHz  (waveguide  designation  WR-4)  to  permit 
measurement  and  characterization  using  commercially  available  instrumentation.  Furthermore, 
the  circuit  is  designed  to  yield  a  180  phase  shift  over  its  operating  band,  a  requirement  for 
efficient  sideband  generation. 

The  phase-shifter  circuit  consists  of  four  basic  components,  i.e.,  (1)  a  pair  of  waveguide-to- 
microstrip  transitions;  (2)  a  hybrid  coupler;  (3)  a  pair  of  varactor  diodes;  and  (4)  a  pair  of 
impedance  transformers.  Fig.  16(a)  shows  a  photograph  of  the  completed  phase-shifter  circuit 
mounted  in  a  WR-4  (170-260  GHz)  waveguide  housing.  Dimensions  (in  micrometers)  of  the 
various  components  are  given  in  Fig.  16(b).  The  system  characteristic  impedance  is  100  Ohms 
and  the  impedance  transformer  between  the  diodes  and  hybrid  is  designed  to  have  an  impedance 
ratio  of  2.  Fused  quartz  is  chosen  as  the  substrate  material  due  to  its  low  cost,  rigidity,  relatively 
low  dielectric  constant,  and  to  permit  monolithic  integration  of  GaAs  varactor  diodes  using  a 
substrate  bonding  process  developed  at  the  University  of  Virginia. 

Performance  of  the  integrated  phase  shifter  was  characterized  and  evaluated  using  an  HP8510C 
vector  network  analyzer  with  Oleson  Microwave  Faboratories  (OMF)  Inc.’s  millimeter-wave 
extension  units.2  The  OMF  extension  units  cover  the  140-325-GHz  range  in  two  standard 
waveguide  bands,  WR-5  (140-220  GHz)  and  WR-3  (220-325  GHz).  Consequently,  the 
measurements  presented  below  represent  a  concatenation  of  the  measurements  performed  over 
both  bands.  It  should  be  noted  that  these  measurements  include  the  small  junction  mismatch 
between  the  WR-3  and  WR-5  flanges  of  the  OMF  extension  modules  and  the  WR-4  flange  of  the 
phase  shifter.  Based  on  simulations  using  HFSS,  the  return  loss  associated  with  this  mismatch  is 
estimated  to  be  no  worse  than  20  dB. 
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Figure  16.  (a)  Photograph  of  the  phase-shifter  chip  mounted  to  the  waveguide  hlock.  (h)  Diagram  of 
the  phase-shifter  chip  layout. 

Fig.  1 7  compares  the  measured  response  of  the  circuit  with  the  smallest  phase  variation  with  that 
predicted  using  Agilent’s  ADS  software.  Over  a  range  of  195-250  GFlz,  the  measured  phase  error 
(referred  to  180  degrees)  is  less  than  15°.  Differences  between  the  measured  and  predicted  phase 
response  for  the  circuit  can  be  attributed  to  several  factors,  including  variations  between  the 
Schottky  diodes  and  errors  in  placement  of  the  circuit  module  within  the  waveguide  channel.  Fig. 
18  shows  the  insertion  and  return  losses  of  the  phase  shifter  measured  from  190  to  250  GFIz.  The 
average  insertion  loss  is  approximately  5.5  dB  over  a  40-GFlz  band  (200-240  GFlz)  and  the 
amplitude  imbalance  between  the  “on”  and  “off’  states  is  better  than  0.5  dB  between  220-240 
GHz. 
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Figure  17.  Measured  phase  shift  between  the  diode  “on”  and  “off’  states  from  190  GHz  to  250  GHz. 
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